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IHIBODOCTIOM 


During  the  operation  of  any  piece  of  mechanical  equip¬ 
ment,  a  certain  amount  of  abrasion  occurs  between  metal 
parts  ccming  in  contact  with  one  another.  The  result  of 
this  abrasion  is  the  deposition  of  minute  quantities  of 
metal  (wearmetals)  in  the  lubricating  oil  of  the  equipment. 

It  was  recognized  by  the  railroais  as  far  back  as  the 
1940's  that  monitoring  the  concentration  cf  wearmetais  in 
the  lubricating  oil  could  be  useful  in  determining  the 
mechanical  condition  of  diesel  engines  and  forewarn  the 
breakdown  of  engine  components. 

Later  in  the  mid  1950' s,  the  Navy  began  a  trial  moni¬ 
toring  pregram  at  the  Naval  Air  Rework  Facility  in  Pensacola 
to  determine  if  oil  analysis  techniques  could  be  useful  for 
detecting  the  abnormal  operation  of  aircraft  engines. 
Following  the  success  of  the  Navy  effort,  the  Army  first 
applied  the  techniques  to  the  monitoring  of  reciprocating 
engine  helicopters  in  1959.  An  evaluation  of  oil  analysis 
programs  was  initiated  by  the  Air  Force  in  1962  and  resulted 
in  the  establishment  of  an  Air  Force  oil  analysis  testing 
program  in  1964. 

Consolidation  of  Department  of  Defense  oil  analysis 
programs  was  initiated  in  1969  with  the  DOD  Equipment  Oil 
Analysis  Frcgram  (ECAP)  and  later  the  Joint  Oil  Analysis 
Program  (JOAP)  .  The  purposes  of  this  consolidation  were  to 
effect  uniformity  and  standardization  in  oil  analysis 
program  equipment,  standards  and  techniques,  consolidation 
cf  procurement  requirements  for  oil  analysis  equipment  and 
the  centralization  cf  responsibility  for  technical  manage¬ 
ment  ana  overall  program  surveillance. 
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The  basic  assumptions  of  tha  spectrometric  analysis  of 
wearmetals  is  that  during  the  normal  operation  of  a  piece  of 
equipment,  the  interacting  of  parts  results  in  the  produc¬ 
tion  of  minute  quantities  of  metals.  This  wearmetal 
production  cccurs  at  a  constant  rate  (as  a  function  cf  time) 
for  properly  operating  pieces  of  equipment.  However,  as  a 
part  enters  into  its  failure  mode,  a  much  higher  level  of 
wearmetal  particles  may  be  deposited  into  the  lubricating 
oil.  An  illustration  of  this  process  for  a  hypothetical 
piece  of  equipment  is  shown  in  Figure  1.1. 

i - 1 
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Figure  1.1  Hypothetical  Wearmetal  Production. 


There  are  limitations  tc  the  types  of  failure  which  can 
be  identified  through  spect rcmetric  analysis.  The  types  cf 
failure  which  can  be  detected  are  failures  which  are  accom¬ 
panied  by  wearmetal  production  and  which  occur  at  a  slow 
enough  rate  to  allow  maintenance  action  after  identification 
of  a  possible  failure.  Some  of  the  types  of  unidentifiable 
failure  are: 


1)  failures  which  produce  large  wearmetal  particles  visi¬ 
ble  tc  the  naked  eye 


2)  failures  which  occur  too  rapidly  to  be  detected  by  the 
current  oil  analysis  procedures 

3)  failures  which  occur  without  appreciable  production  of 
wearmetals 

To  identify  impending  failure  conditions,  wearmetal 
levels  in  an  oil  sample  from  the  eguipment  studied  are 
determined  through  spect  rom  etric  analysis.  Next,  the  labora¬ 
tory  personnel  calculate  a  10  hour  PPM  increase  figure  by 
finding  the  change  in  PPM  between  samples,  dividing  that 
figure  by  the  hours  between  samples  (giving  a  rate  of 
increase  per  hour)  and  multiplying  by  10  (giving  the  rate  of 
increase  for  10  hours).  A  guideline  chart  giving  wearmetal 
PPM  criteria  for  normal,  marginal,  high  and  abnormal  ranges 
and  for  abnormal  10  hour  trends  is  consulted  for  the  appli¬ 
cable  type  of  equipment.  These  guidelines  are  not 
recommended  wearmetal  values,  but  indicate  ranges  of  values 
within  which  the  PPM  levels  of  similar  equipment  have 
historically  fallen. 

A  determination  is  made  as  to  whether  the  amount  of 
wearmetals  in  the  sample  is  considered  in  the  normal, 
marginal,  high,  or  abnormal  range  and  whether  the  PPM 
increase  in  the  10  hour  period  is  abnormal.  Based  on  the 
range  of  the  sample,  the  range  of  the  previous  sample,  the 
trend  reading  and  whether  the  oil  sample  was  a  routine 
submission  or  the  result  of  a  laboratory  request,  a  recom¬ 
mendation  is  made  utilizing  the  Decision  Making  Guidance 
Table  (Table  6-1)  of  TM  38-301.  Basically  the  table  calls 
for  increased  sampling  if  there  is  an  indication  of  abnormal 
operations.  Continued  abnormal  readings  result  in  the  dead¬ 
lining  of  the  equipment  and  a  recommendation  that 
maintenance  action  be  initiated  to  determine  the  cause  of 
the  high  level  of  wearmetals.  In  some  cases  a  combination 
of  abncrmal  wearmetal  readings  can  be  used  to  pinpoint  the 


specific  area  of  the  equipment  which  is  probably  in  the 
process  of  failing. 

The  maintenance  personnel  must  than  provide  followup  to 
the  oil  analysis  laboratory  stating  whether  rheir  indication 
of  impending  failure  was  correct  (known  as  a  'hit')  or 
incorrect  (known  as  a  'miss*)« 

Thus,  through  the  combined  efforts  of  the  oil  analysis 
laboratory  and  the  maintenance  personnel,  many  impending 
failures  can  be  detected  and  corrected,  preventing  further 
damage  (cr  perhaps  destruction)  of  the  equipment  and  serious 
injury  to  the  personnel  operating  ir. 


II.  NAT ORE  OF  THE  PROBLEM 


Ext5nsi.cn  of  the  oil  analysis  program  to  Army  ncn- 
aeronautical  equipment  is  fairly  recent,  therefore  there  are 
a  number  of  engines  in  the  Army  inventory  which  are  still 
lacking  wearmetal  and  trend  guidelines.  The  DARCOM  Materiel 
Readiness  support  Activity  at  Lexington,  Ky.  provided  a  list 
of  engines,  used  in  Army  ground  equipment,  which  were 
lacking  published  guidelines.  Of  the  engines  listed,  the 
one  which  was  most  widely  used  (in  terms  of  types  of  equip¬ 
ment  of  different  nomenclature)  was  the  Continental 
LD/LDS/LCT  465/465-1  (henceforth  referred  to  as  the 
IDS- 465 ) .  This  is  a  six  cylinder  diesel  engine  which  powers 
a  large  number  of  the  Army’s  2  1/2  Ton  and  5  Ton  trucks.  A 
list  of  the  types  of  vehicles  powered  by  the  engine  is  at 
Appendix  A. 

The  goal  of  the  analysis  was  to  examine  historical  data 
(from  DD  Form  2027-  Cil  Analysis  Record)  on  a  large  sample 
of  different  vehicles  operated  in  different  parts  of  the 
country,  which  were  equipped  with  the  LDS-465  engine.  From 
this  analysis  it  was  hoped  to  derive  guidelines  for  what 
could  be  considered  normal,  marginal,  high  and  abnormal 
concentrations  of  wearmetals  in  the  lubricating  oil.  These 
guidelines  were  to  be  general  (i.e.  applicable  to  all  types 
cf  equipment  using  the  subject  engine)  in  keeping  with  the 
format  of  the  guidelines  published  in  TM  38-301  and  TM 
38-301-1.  in  addition,  it  was  desired  to  determine  the 
relationship  between  operating  hours  of  the  piece  of  equip¬ 
ment  and  the  rate  of  wearmetal  production.  This  information 
was  to  be  used  to  develop  guidelines  on  what  was  to  be 
considered  an  abnormal  increase  in  wearmetal  PPM  over  a  10 
hour  operating  period. 

1 1 
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Due  tc  time  constraints  in  conducting  the  analysis,  Army 
wide  data  could  not  be  obtained.  The  scope  of  the  analysis 
was  therefore  limited  to  a  sample  of  data  on  equipment 
serviced  by  the  Ft.  Ord  Spectrometric  Oil  Analysis 
Laboratory.  This  limitation  on  the  scope  of  the  analysis 
also  placed  a  limitation  on  the  applicability  of  the  results 
in  that  they  are  only  representative  of  the  equipment  at  Ft. 
Ord  and  are  based  on  a  relatively  small  sample,  thus  exten¬ 
sion  of  the  results  to  Army  wide  equipment  should  be 
avoided.  The  methodology,  however,  may  be  useful  in 
analyzing  a  much  larger  set  of  data  in  order  to  provide 
guidelines  applicable  to  Army  wide  equipment. 
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III.  AN ALT  SIS  METHODOLOGY 


A.  DATA  COLLECTION 

Historical  oil  analysis  records  (DD  Form  2027)  for  the 
period  January  1980  through  January  1983  were  examined  for 
76  vehicles  from  Active  Army  units  located  at  Ft.  Ord.  The 
records  studied  were  for  2  1/2  Ton  Cargo  trucks  of  the  type 
M35A2  or  M35A2C.  A  list  of  serial  numbers  for  the  vehicles 
examined  is  at  Appendix  Q. 

The  M35A2  and  M35A2C  are  tactical  vehicles  and  as  such, 
the  guidance  on  routine  sampling  intervals  is  once  every  two 
months.  Samples  are  taken  while  the  oil  is  warm  either  by 
the  use  cf  a  dip  tube  or  by  a  drain  outlet.  Detailed 
descriptions  of  the  sampling  methods  are  found  in  TM  38-301 
cr  TB  43-0211. 

An  Oil  Analysis  Bequest  (DD  Form  2026)  is  filled  cut 
giving  such  information  as  equipment  identification  data, 
cime  cil  sample  taken  and  method  used  (drain  or  dip  tube)  , 
hours  since  oil  change,  hours  since  overhaul,  and  any  state¬ 
ments  abcut  unusual  operation  of  the  equipment.  The  samples 
are  then  sent  to  the  Oil  Analysis  Laboratory  at  Ft.  Ord 
where  a  spectrometric  analysis  is  performed  to  determine 
wearmetal  levels  in  the  sample. 

The  analysis  is  accomplished  by  the  use  of  the  Baird 
Emission  Spectrometer  Model  A/E350-3.  A  small  portion  of 
the  oil  sample  is  placed  in  the  spectrometer.  A  carbon  disk 
electrcde  rotates  through  the  sample,  picking  up  a  film  of 
oil  which  is  burned  by  an  electric  arc  between  a  fixed  elec¬ 
trode  and  the  rotating  disk.  The  energy  radiated  by  the 
burning  oil  is  separated  into  its  component  wavelengths  and 
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the  intensity  measured.  This  is  compared  to  the  intensity 
measured  from  calibration  samples  to  determine  the  parts  per 
million  (PPM)  levels  of  the  various  wearmetals  in  the 
sample.  There  is  a  degree  of  variablility  associated  with 
the  spectrometer  reading.  TM  38-301  states  the  spectrometer 
has  a  tolerance  of  ±1.0  PPM;  however  it  seems  more  likely 
the  variability  is  a  function  of  the  magnitude  of  the 
reading  as  suggested  by  the  spectrometer  manufacturer. 

After  spect rometric  analysis,  the  following  information 
is  then  entered  into  the  DD  Form  2027  (Oil  Analysis  Record) 
for  the  sampled  piece  of  equipment. 

1.  Sample  number  (laboratory  assigned) 

2.  Data  index  (  data  correction  or  maintenance  feedback) 

3.  Julian  date  sample  analyzed 

4.  Laboratory  response  time 

5.  Last  laboratory  recommendation 

6.  Hcurs  since  overhaul 

7.  Hcurs  since  oil  change 

8.  Reason  for  the  sample 

9.  Wearmetal  levels 

10.  Post  analysis  data  (  lab  recommendations  and  mainte¬ 
nance  followup) 

None  of  the  records  studied  indicated  laboratory  recom¬ 
mendations  for  maintenance  action  and  thus  no  reports  of 
maintenance  followup  existed.  It  was  also  found  that 
approximately  25*  of  the  historical  records  lacked  data  for 
hours  since  oil  change  and  hours  since  overhaul.  This 
precluded  any  systematic  analysis  of  the  entire  sample  of 


wearmetal  levels  as  a  function  of  operating  hours.  As  a 
result,  to  determine  wearmetal  guidelines,  an  analysis  was 
conducted  which  did  net  require  operating  hour  information. 

The  suggested  oil  sampling  intervals  were  also  not 
adhered  to  in  all  cases.  As  a  result,  the  number  of  records 
available  for  each  equipment  during  the  time  period  studied, 
ranged  from  14  to  17  samples. 

B.  DETER HI NATION  OF  8BARMBTAL  LEVEL  GUIDELINES 

Faced  with  these  constraints,  the  method  of  analysis 
chosen  was  to  look  at  all  of  the  PPH  levels  available  for 
the  76  vehicles,  plot  a  frequency  distribution  and  examine 
the  quantiles  from  the  empirical  distribution. 

The  director  of  the  Ft.  Ord  Oil  Laboratory  provided 
guidance  on  the  major  wearmetals  of  concern  for  the  LDS-465 
engine.  They  were: 

1.  Iren  <Fe) 

2.  Aluminum  (Al) 

3.  Chromium  (Cr) 

4.  Copper  (Cu) 

5.  Lead  (Pb) 

Thus  the  scope  of  the  study  was  limited  to  these  five 
wearmetals.  The  wearmetal  data  was  read  by  the  FORTRAN 
program  WEARMTL  (see  Appendix  F  for  computer  program  list¬ 
ings)  which  produced  an  output  frequency  histogram  for  each 
wearmetal  along  with  a  number  of  sample  statistics.  In  all, 
a  total  of  1255  sample  points  were  available.  The  60th, 
70th,  80th,  90th,  95th,  96th,  97th,  98th,  and  99th  quantiles 
were  determined  from  the  data.  No  decisions  were  made  on 


what  labels  (i.e.  normal,  marginal,  high,  or  abnormal)  to 
place  on  these  quantiles,  but  rather  to  provide  this  infor¬ 
mation  to  the  Oil  Analysis  Laboratory  director  for  use  in 
his  own  subjective  judgement.  The  quantiles  found  are 
listed  in  Table  I.  Output  histograms  are  at  Appendix  C. 


In  addition  to  wearmetal  levels,  a  major  concern  in  the 
oil  analysis  program  is  wearmetal  trends.  Trend  guidelines 
are  established  for  a  standard  10  hour  period  and  give  what 
is  considered  to  be  an  abnormal  increase  in  wearmetal  levels 
for  that  period.  In  the  study  performed  by  the  ARINC 
Corporation  for  the  0.  S.  Air  Force  [Ref.  1,  pp  10-14]  the 
method  of  piecewise  or  segmented  regression  was  used  to 
determine  the  point  of  onset  of  failure.  In  this  method, 
wearmetal  concentration  data  on  a  number  of  similar  pieces 
of  equipment  is  divided  into  two  groups,  that  data  repre¬ 
senting  normal  wear  and  that  representing  abnormal  wear. 
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The  hours  associated  with  each  sample  are  converted  from 
hours  since  oil  change  to  hours  prior  to  detection  of 
failure.  This  indicates  the  necessity  of  maintenance  feed¬ 
back  to  determine  or  verify  a  point  of  failure  detection. 
Straight  lines  are  then  fit  through  both  sets  of  data  using 
the  procedure  detailed  by  Hudson  [Ref.  2].  The  point  of 
intersection  of  the  two  lines  is  taken  to  be  the  tine  prior 
to  detection  at  which  failure  commenced  (see  Figure  3.1). 
The  slope  of  the  line  following  this  point  could  then  be 
used  to  calculate  abnormal  10  hour  trends. 


i -  - - 1 


Figure  3. 1  Fitting  Data  By  Piecewise  Regression. 

Since  no  maintenance  feedback  was  available,  a  technique 
was  proposed  which  would  indicate  abnormally  high  metal 
trends  without  requiring  maintenance  information  on  actual 
failure  conditions.  This  method  assumes  that: 

1.  Wearmetal  levels  may  depend  on  operating  hours  of  the 
equipment. 

2.  A  possibly  different  relationship  exists  for  the  nor¬ 
mal  and  abnormal  periods  of  operation. 
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3.  The  rata  of  wearmetal  production  is  higher  during 
failure  than  during  normal  operation. 

Instead  of  fitting  piecewise  regression  lines  to  the 
data,  the  proposed  methodology  treats  the  slope  of  each  line 
segment  (whose  end  points  are  formed  by  successive  wearmetal 
readings)  as  an  estimate  cf  the  rata  of  wearmetal  produc¬ 
tion.  An  illustration  of  this  tachnigue  is  given  in  Figure 
3.2.  Each  dotted  line  segment  provides  an  estimate  of  the 


Figure  3.2  Piecewise  Estimation  10  Hour  Wearmetal  Trends. 

production  rate.  If  the  wearmatal  production  rate  is 
linear,  the  relationship  is  of  tha  form: 

Y(i)  =  B  (0)  ♦  B  ( 1)  *X  (i)  +  5(i) 


where : 


X  (i)  3  hours  since  oil  change  for  the  ith  sample 


I  (i)  =  wearmetal  reading  for  the  ith  sample 
B  (0)  =  intercept  of  the  regression  line 


E  (i)  =  deviation  of  the  actual  wearmetal  reading 
from  the  regression  line. 

The  above  definitions  apply  throughout  the  rest  of  the 
thesis. 

The  slcpe  of  each  successive  line  segment  is  found  by  the 
formula 

slope  =  (T  (i+1)  -T (i) )  /(X  (i+ 1)  -X  (i) ) 
substituting  the  regression  equations: 

Slope  =  B  (0)  +  (B  (1)  *X  (i+1))  +  E  (i  +  1)  -B  (0)  -  (B(1)*X(i)  )-E  (i) 

X(i  +  1)  -X(i) 

=  £(  1)  ♦  (E  (i  +  1 )  -E(i) )  /  (X  (i+1)  -X  (i) ) 

Since  the  regression  slope  B(1)  is  constant  for  any  given 
set  cf  data,  each  pairwise  slops  will  deviate  from  the 
actual  slcpe  by  a  factor 

(E  (i+1)  -E(i) )  /  (X  (i+1)  -X  (i)  ) 

Thus,  the  'goodness'  of  this  method  is  determined  by  the 
magnitude  of  the  variance  in  the  residuals  and  the  length  of 
time  between  wearmetal  samples.  Samples  wizh  small  vari¬ 
ances  in  their  error  terms  and  large  sampling  intervals  will 
approximate  the  regression  slope  much  better  than  those  with 
large  variances  and  small  sampling  intervals. 

The  next  step  in  the  methodology  is  to  convert  each  of 
these  slcFes  to  a  PPB  increase  for  a  10  hour  period.  This 
is  done  by  multiplying  the  slope  by  10: 

Increase  in  PPB  =  1 0 *  (B ( 1)  ♦  ( ( (E  (i  +  1)  -  E  (i) )  /  (X  (i+  1 ) -X  (i) ) )  ) ) 
in  10  hours 

Thus  each  of  the  slope  segments  (regardless  of  the  time 
interval  between  the  samples  )  is  converted  to  a  change  in 
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PPM  which  would  be  expected  during  any  10  hour  period  within 
this  interval.  The  10  hour  PPM  rate  of  change  measurements 
are  then  grouped  together  in  a  frequency  histogram.  If  the 
rate  of  wearmetal  production  is  higher  during  failure  than 
during  normal  operation,  then  the  upper  quantiles  of  the  EPM 
rate  of  change  distribution  should  represent  these  samples 
from  the  failure  mode  of  the  equipment. 

This  methodology  was  applied  to  the  Fort  Ord  data  (for 
those  data  elements  which  had  hours  since  oil  change  infor¬ 
mation  available)  through  the  use  of  the  FORTRAN  program 
ANALYSIS  (Appendix  F) .  The  logic  used  in  computing  the 
slopes  was: 

1.  If  no  hours  since  oil  change  data  was  available,  the 
sample  was  discarded. 

2.  If  hours  since  oil  change  was  greater  than  the  pre¬ 
vious  reading,  the  slope  was  determined  by: 

(Y  (i  +  1)-Y(i) )  / (X  (i+1)  -X  (i)  ) 

3.  If  hours  since  oil  change  was  less  than  that  of  the 
previous  sample,  it  was  assumed  an  oil  change  took 
place  and  thus  the  slope  was  calculated  by: 

Y(i)/X(i) 


Since  approximately  25X  of  the  records  lacked  time  since 
oil  change  data,  only  8  76  tread  values  were  computed. 
Histograms  of  the  trend  frequencies  are  at  Appendix  D.  The 
results  obtained  showed  a  large  variance  in  the  PPM  rate  of 
change  levels.  This  may  have  been  due  to  errors  in  data 
reporting  since  a  large  factor  in  the  variance  was  the  exis- 
tance  of  a  few  extremely  large  or  small  measurements.  As  in 
the  case  of  the  wearmetal  levels,  selected  quantiles  were 
obtained  and  are  shown  in  Table  II. 


TABLE  II 

Selected  10  Hoar  Tread  Quantiles 


Wearmetal 

Iren 

Aluminum 

Copper 

Chromium 

Lead 


Quantile 


6  0th 

70th 

80  th 

90xh 

95th 

6.0 

9.3 

15.7 

31.9 

66.9 

1.2 

2.0 

3.5 

8.5 

20.0 

1.1 

1.9 

3.3 

7.5 

15.9 

0.5 

0.9 

1.5 

3.3 

8.0 

2.4 

4.5 

8.8 

18.  3 

37.5 

96th  97xh  98th  99th 

86.3  140.  190.  230. 

23.8  30.0  41.3  70.0 

20.9  26.9  35.0  45.9 
10.0  13.0  20.0  23.3 

42.4  52.5  80.0  110. 


IV.  IlIOSTBATIOM  OP  THE HD  ANALYSIS  TECHHIOOES  OSING 

SIHOIATBD  DA JA 

The  technique  proposed  for  analyzing  trends  was  not 
adequately  illustrated  by  the  actual  data  for  the  76  vehi¬ 
cles  in  the  sample  due  to  the  suspect  validity  of  the  hours 
since  oil  change.  To  further  illustrate  the  technique,  a 
set  cf  data  was  generated  which  fit  the  assumptions  set 
forth  when  proposing  this  technique.  Those  assumptions 
were: 

1.  PFH  levels  may  be  related  to  the  operating  hours 
since  oil  change. 

2.  The  wearmetal  production  process  may  be  different  for 
the  normal  and  abnormal  modes  of  operation. 

3.  The  rate  of  wearmetal  production  is  higher  during 
failure  than  during  normal  operation 

Thus  simulated  wearmetal  data  was  generated  by  assuming: 

Y  (i)  =  B  (0)  +  B(1)*X(i)  ♦  E  (i) 

Where  the  variables  are  defined  as  in  chapter  III. 

Since  the  proposed  technique  approximating  the  regression 
line  slope  by  the  individual  piecewise  slopes  was  indepen¬ 
dent  of  the  intercept  term  (i.e.  B(0)),  it  was  only 

necessary  to  derive  estimates  for  B(1)  and  the  variance  of 
the  error  terms. 

To  determine  a  representative  value  to  use  for  B(1), 
records  fcr  the  76  trucks  were  examined  to  determine  whether 
any  of  the  vehicles  had  wearmetal  data  for  iron  which  exhib¬ 
ited  the  following  characteristics: 


1.  Continuous,  relatively  uniform  increases  in  hours 
since  oil  change  data 

2.  No  apparent  gaps  in  the  data 

3.  No  abnormal  drops  in  EPM  levels  without  an  accompa¬ 
nying  record  of  oil  change 

4.  At  least  8  data  points  fitting  the  above  requirements 


A  total  of  ten  vehicles  met  the  above  criteria.  The  data 
from  these  vehicles  was  entered  into  the  APL  program  HEGPESS 
[aef.  3],  which  fit  a  least  squares  linear  regression  line 
to  the  data.  A  plot  of  the  residuals  was  examined  to  deter¬ 
mine  if  there  were  any  apparent  violations  of  the  assumption 
of  normality.  In  addition,  the  fi2  values  were  examined  to 
determine  if  there  was  a  strong  relationship  between  the 
dependent  (PPM  levels)  and  independent  (hours  since  oil 
change)  variables.  Those  samples  which  had  32  values  of  .9 
or  greater  and  did  not  violate  the  subjective  test  for 
normality  of  residuals  were  chosen.  A  total  of  four  vehi¬ 
cles  met  these  criteria.  The  iron  3(1)  values  for  these 
vehicles  are  shown  in  Table  III. 


TABLE  III 

Regression  Line  Slopes 


ehicle  number 

B(1) 

3903458 

.  364  8 

3818363 

.  4277 

3 85 C 984 

.514  1 

3903880 

.3306 
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This  technique  is  not  meant  to  represent  a  statistically 
valid  method  for  determining  the  behavior  of  the  wearmetal 
data.  It  is  meant  only  to  provide  a  general  idea  of  the 
magnitude  of  the  wearmetal  production  rates. 

The  variance  was  assumed  to  be  1  for  the  sake  of 
simplicity. 

Thus  to  generate  the  data,  all  that  was  necessary  was  to 
determine  the  time  at  which  the  sample  was  to  be  taken, 
multiply  it  by  the  slope  ( B  ( 1 ) )  and  add  on  an  error  term 
wnich  was  distributed  N  (0,  1)  .  For  illustrative  purposes, 
times  were  generated  at  exact  ten  hour  intervals  (i.e.  10, 

20,  30  etc.  hours)  .  The  error  terms  were  generated  using 

the  IIS  A  NEON  II  [Bef.  4]  random  number  generating  package. 
To  simulate  abnormal  operations,  data  was  generated  which 
had  a  slcpe  twice  that  of  the  normally  operating  equipment. 
Error  terms  were  generated  in  the  same  manner  as  for  the 
normally  operating  equipment.  The  data  generated  for  the 
four  pieces  of  equipment  were  combined  and  a  frequency 
distribution  (histogram)  was  formed  through  use  of  the 
FORTRAN  program  ANAIGEN  (Appendix  F)  .  Five  hundred  fifty 
•normal'  data  elements  and  fifty  'abnormal*  data  elements 
were  generated  for  each  of  the  four  values  of  B(l).  The 
cutput  frequency  histogram  is  at  Appendix  E. 

Upper  quantiles  cf  the  distribution  were  determined  as 
was  dene  with  the  wearmetal  levels.  The  quantiles  found  are 
in  Table  IV. 

A  problem  may  arise  any  time  data  from  different  equip¬ 
ment  is  grouped  together.  Unless  the  assumption  that  all 
equipment  of  the  same  type  operate  in  the  same  manner  is 
true,  one  could  identify  a  measurement  as  abnormal  in  terms 
cf  the  entire  population  when  in  fact  it  is  within  the 
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TABLE  IV 

Selected  Quantiles  of  10  Hour  tfearaetal  Production 


Q  OANTILE 

WEARMETAL  60tn  70th  80th  90th  95th  96th  97th  98th  99th 
Iren  4.69  5.22  5.84  6.97  8.1  8  3.83  9.  46  10.14  10.  81 


normal  operating  range  of  the  equipment  from  which  the  meas¬ 
urement  was  taken.  To  preclude  this  happening  it  would  be 
preferable  if  guidelines  could  be  generated  for  each  equip¬ 
ment  based  cn  the  operating  history  of  that  equipment  alone. 

An  example  of  the  benefits  of  this  method  is  shown  using  the 
generated  data  from  vehicle  3SQ3458. 

Since  the  error  terms  generated  ware  distributed  Normal 
(0,1),  the  distribution  of  the  PPM  reading  for  each  gener-  fc 

atea  sample  should  have  been  Normal  (10*B(1),2).  This 
follows  since  the  generated  intervals  were  10  hours  there¬ 
fore  each  10  hour  PPM  reading  would  be  10  times  the  slope  of 
the  regression  line  (which  was  constant)  plus  the  difference 
in  the  pairwise  consecutive  error  terms.  Since  the  mean  of 
the  error  terms  was  0,  the  mean  of  the  PPM  readings  would  be 
0+10*E(1),  while  the  variance  would  be  -he  variance  of  B(1)r 
which  was  zero,  plus  the  variance  of  the  difference  in  the  fcj 

error  terms  which  was  2.  Thus  it  would  be  expected  that  the 
frequency  distribution  for  each  generated  sample  would  be  a 
combination  of  two  normal  distributions  offset  from  each 
other  by  10  times  the  difference  in  their  slopes  (means).  A  £ 

histogram  of  the  data  generated  from  the  B(1)  value  for 
vehicle  3803458  (Appendix  E)  shows  this  to  be  the  case  by 
the  distinct  bimodality  of  the  distribution. 

» 

25 


Thus  if  tha  assumptions  of  the  vearmetal  bahavicr  during 
normal  and  abnormal  operations  hold,  and  the  difference  in 
wearmetal  production  rates  during  normal  and  abnormal  opera¬ 
tions  is  large  in  comparison  to  tha  variance  in  the  sample 
cf  normal  readings,  then  it  should  be  quite  apparent  whan  a 
reading  falls  out  of  the  range  of  normal  variability. 
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7.  CQNCIOSIONS  AND  BECOHHEMDATIONS 

The  determination  of  quantiles  from  the  empirical 
distributions  of  wearmetal  levels  and  of  wearmetal  trends 
was  done  to  provide  the  director  of  the  Ft.  Ord 
Spectrcmetric  Oil  Analysis  laboratory  with  general  guide¬ 
lines  on  what  ranges  of  values  he  could  expect  for  vearmetal 
levels  and  trends  for  the  Continental  LD/LDS/LDT  465/465-1 
based  on  historical  data  for  equipment  supported  by  his 
laboratory.  The  presentation  of  results  was  different  than 
that  in  TM  38-301-1  in  that  no  attempt  was  made  to  attach 
labels  (i.e.  normal,  marginal  etc.)  to  these  quantiles. 
This  approach  was  taken  to  provide  more  flexibility  and 
information  with  which  to  make  subjective  judgements  on  the 
condition  cf  a  piece  of  equipment  since  ultimately  any 
recommendation  made  by  an  oil  analysis  laboratory  is  subjec¬ 
tive  based  not  only  cn  published  guidelines  for  a  type  of 
equipment,  but  also  a  knowledge  of  the  operational  history 
of  the  individual  equipment  in  question. 

The  fact  that  operational  history  is  an  important  factor 
in  the  recommendations  made  by  the  laboratory  indicates  the 
current  practice  of  griping  all  equipment  of  the  same  type 
together  and  applying  general  guidelines  to  the  composite 
sample  may  not  in  fact  be  adequate  for  all  individual  pieces 
of  equipment,  being  too  conservative  for  equipment  which 
normally  produce  wearmetal  at  a  relatively  high  rate  or  too 
liberal  for  equipment  which  produce  wearmetal  at  a  rela¬ 
tively  lew  rate.  This  problem  is  illustrated  in  Figure  5.1. 
In  this  figure,  samples  from  two  hypothetical  pieces  of 
equipment  are  shown.  A  regression  analysis  is  performed  on 
each  sample  individually  giving  regression  lines  1  and  2 


REGRESSION  OF  MULTIPLE  SAMPLES 


TIME  (HOURS) 


respectively.  These  lines  are  constructed  to  minimize  the 
sum  cf  the  squared  deviations  of  each  sample  point  from  the 
line.  Unless  the  two  regression  lines  are  collinear  (ar. 
extremely  unlikely  situation),  a  regression  analysis 
performed  on  the  combined  data  would  result  in  a  regression 
line  (line  3)  which,  while  being  the  best  fit  for  the 
combined  data,  would  have  a  larger  residual  sum  cf  squares 
than  either  of  the  individual  sample  regression  lines  indi¬ 
cating  the  fit  of  the  regression  line  for  the  combined  data 
was  not  as  good  a  fit  as  that  of  the  individual  samples. 
Additionally,  the  combined  regression  line  would  overesti¬ 
mate  the  wearmstal  production  rate  of  engine  2  and 
underestimate  the  wearmetal  production  rate  of  engine  1. 
Thus  a  guideline  based  on  the  combined  samples  might  fail  to 
identify  an  impending  failure  in  engine  2  or  erroneously 
signal  an  impending  failure  in  engine  1. 

This  same  problem  applies  to  the  methodology  proposed  of 
estimating  the  regression  line  by  the  individual  10  hour 
trend  readings.  While  each  10  hour  trend  reading  could  be 
considered  a  fairly  good  estimate  of  the  slope  of  the 
regression  line  for  that  individual  piece  of  equipment,  it 
would  likely  be  a  poor  estimate  of  the  slope  of  the  regres¬ 
sion  line  obtained  by  combining  the  data  from  all  pieces  of 
equipment  of  the  same  model.  Again,  the  wearmetal  produc¬ 
tion  characteristics  of  the  individual  equipment  would  be 
masked  by  the  combined  characteristics  of  the  other  pieces 
cf  equipment. 

while  developing  individual  guidelines  on  each  piece  of 
equipment  is  much  mere  complicated  than  setting  general 
guidelines  for  ail  equipment  of  the  same  type,  these  indi¬ 
vidual  guidelines  would  allow  the  laboratory  director  to 
make  a  much  more  accurate  recommendation  after  analyzing  an 
oil  sample. 


I  up  1 6  id  sn  ration  of  this  methodology  at  the  laboratory 
level  would  require  computational  capability  (i.e.  a 
computer) .  Since  the  distribution  of  the  10  hour  PPM  trends 
for  an  individual  piece  of  equipment  is  expected  to  be 
normal,  any  desired  quantiles  of  the  distribution  could  be 
derived  from  the  sample  mean  and  standard  deviation  and  the 
standard  normal  tables.  As  additional  data  elements  were 
added  (as  a  result  of  additional  samples)  these  statistics 
would  provide  an  increasingly  better  estimation  of  the  true 
wearmetal  trend  characteristics  of  the  individual  equipment. 

The  goodness  of  any  predictive  methodology  is  related  to 
the  quality  of  the  data  on  which  the  prediction  is  made. 
Thus  to  derive  the  maximum  benefits  from  the  oil  analysis 
program,  additional  emphasis  must  be  placed  on  the  proper 
collecting  of  oil  samples  and  reporting  of  equipment  oper¬ 
ating  data.  This  will  require  not  only  increased  command 
emphasis,  but  convincing  maintenance  personnel  of  the 
benefits  of  the  program  in  reduced  maintenance  expense, 
increased  readiness  and  greater  safety. 


ASISSSI1  k 

TIPE  VEHICLES  USING  THE  CONTINENTAL  LD/LDS/ LDT-465  ENGINE 


END  ITEM  NOMENCLATURE 


M34A2 

2  1/2 

M35A2 

•1 

M35A2C 

II 

H36A2 

II 

M45A2 

2  1/2 

M46A2 

2  1/2 

M49A2C 

2  1/2 

M50A2 

2  1/2 

M50A3 

2  1/2 

H51A2 

5  Tor. 

M52A2 

5  Ton 

M54A2 

5  Ton 

M54A2C 

II 

M55A2 

II 

M61A2 

5  Ton 

B63A2 

•i 

M109A3T 

2  1/2 

8135  A3 

In  st 

H246A2 

5  Ton 

Ton  Cargo  Truck 

it 

it 

it 

Ton  Bolster,  truck 
Ton  truck  chassis 
Ton  Tank  truck,  fuel 
Ton  Tank  truck,  H20 
Ton  Tank  truck,  H20 
Dump  truck 
Tractor  truck 
Cargo  truck 

If 

II 

truck  chassis 

ii 

Ton  truck,  tractor,  wrecker 
Repair  truck 
truck,  tractor,  wrecker 


AP£ENDi£  B 

LIST  OF  SEBIiL  NOE  BEES  OF  ENGINES  STUDIED 


2  5  88  0  A 

381 7608 

2811 

12688 

3929002 

31244 

3903358 

392  3271 

3903782 

3923514 

390  1901 

3902143 

3929006 

3920268 

3921134 

3914728 

392  1828 

3923959 

3901793 

390  1793 

13778 

3902619 

20279 

6905 

3765142 

3806178 

382045 

3  9  22  7  9  5 

3903772 

38176592 

3827828 

10305783 

980716 

3891  139 

261  97 

3818461 

3862237 

390  3451 

880819 

391707 

390 1904 

3902362 

3902420 

3903137 

24165 

3924675 

3924712 

3902430 

3090  8 

390  1455 

79168 

3903418 

3922157 

75145 

3  9  23  8  4  57 

391  2945 

3914805 

3925101 

390  3884 

3047329 

3903458 

3939020 

3818363 

3901270 

390  3993 

3900086 

8805274 

4076463 

3925768 

3  650984 

3903880 

3903421 

39032  12 

381 6622 

3929195 

3903350 
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HISTOGRAM  OF  IRON  PPM 


PPM  OF  IRON 


HISTOGRAM  OF  ALUMINUM  PPM 
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PPM  OF  ALUMINUM 


HISTOGRAM  OF  CHROMIUM  PPM 
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PPM  OF  CHROMIUM 


HISTOGRAM  OF  COPPER  PPM 


OF  COPPER 


HISTOGRAM  OF  LEAD  PPM 


PPM  OF  LEAD 


APPENDIX  D 

OUTPUT  POS  THEID  ANALYSIS 


The  results  of  this  portion  of  the  analysis  ware  highly 
variable.  The  existance  of  a  few  very  large  and  very  small 
values  prevented  the  display  of  the  bulk  of  rhe  results  ir. 
detail.  To  expand  on  the  display  of  the  non-extreme 

results,  the  10  largest  and  10  smallest  values  were  trun¬ 
cated  prior  to  producing  the  histograms.  The  values 
truncated  for  each  wearmetal  are  listed  in  Table  v. 


TABLE  V 

Truncated  Bearmetal  Trend  Values 


IECN  ALUMINUM  CHROMIUM 

LOW  HIGH  LOW  HIGH  LOW  H 


LOW 

HIGH 

LOW 

HIGH 

LOW 

HIGH 

-680 

223. 33 

-280 

64 

-70 

22.22 

-450 

230 

-190 

70 

-20 

23.33 

-370 

230 

-60 

70 

-12.31 

27.5 

-135.71 

232 

-40 

70 

-10 

30 

-125. 83 

25  0 

-30 

75.71 

-10 

40 

-  100 

263. 33 

-3  0 

110 

-10 

40 

-100 

27  0 

-28.57 

110 

-10 

50 

-92. 5 

383.  33 

-25 

130 

-10 

50 

-86.67 

610 

-23.3  3 

130 

-10 

53.33 

-80 

915 

-22.5 

355 

-6.67 

70 

C  CEPES 

LEAD 

LOW 

HIGH 

LOW 

HIGH 

-100 

45.46 

-380 

103.33 

-50 

46 

-250 

110 

-30 

50 

-251. 43 

1 15 

-21. 43 

62  .5 

-120 

117 

-  19. 23 

70 

-90 

126.67 

-18.33 

80 

-83.33 

130 

-  12. 67 

83.33 

-8  0 

250 

-12.  5 

110 

-60 

252.86 

-12.  5 

160 

-56.67 

275 

-11.52 

180 

-40 

380 
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Selected  statistics  on  the  10  hour  trends  are  shewn  in  Table 

VI. 


TABLE  VI 
Trend  Statistics 


WEARMETAL 


STATISTIC 

Iron 

Alum 

Chrom 

CODD 

Lead 

MEANS 

1  1.88 

3.08 

2.97 

1.32 

5.60 

STE. DEV 

62.69 

20.81 

12.86 

6.1  1 

32.  54 

L.  QRTIS 

0.13 

0.00 

0.00 

0.00 

0.00 

MEDIANS 

3.49 

0.71 

0.61 

0.23 

1.  37 

U.  QRTIS 

1  1 .76 

2.58 

2.44 

1.11 

6.19 
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O  Q 

< 

LlJ 


cr 
3 
o  O 
X 


A0N3n03dJ 
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oc 

x 
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2 
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X 

9 
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10  HOUR  READINGS  FOR  CHROMIUM  (PPM) 
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Kznnn  s 

OUTPUT  FOB  TBEND  ANALYSIS  USING  SIMULATED  DATA 


This  appendix  contains  histograms  of  the  wearmeral  tren 
analysis  for  the  simulated  wearmetal  data.  The  first  histo 
gram  is  for  the  combined  data  of  four  engines  while  th 
second  is  a  histogram  for  an  individual  engine. 
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APPENDIX  F 
COMPUTER  PE 0 GRABS 


C  ****  EBGGBAM  WEARMTL**** 

C  PROGRAM  READS  IN  WEARMETAL  LEVELS  FOR  DIFFERENT 
C  EQUIPMENT  AND  PLOTS  FREQUENCY  HISTOGRAMS 
C  RICHABE  F.  BAUER  11  APR  1983 

REAL  HR  (126  0)  ,FE(126  0)  ,  AL  (1260)  ,CR(1260)  f  CU  (1260) 
♦  PEM260)  ,  INDX  (1260) 

INTEGER  I,J 
DIMENSION  NAME  (200) 

J  =  C 
1*0 

DO  15  J=1  ,75 

READ  (5,100) 

WRITE  (6  ,270) 

READ  IN  THE  WEARMETAL 
10  CONTINUE 
1=1+1 

READ  (5,110)  HR  (I)  ,FE  (I)  ,AL(I)  ,C3(I)  ,CU  (I),  PB  (I) 
WRITE  (6, 250)  HR  (I)  ,  F  E  (I)  ,  AL  (I )  ,CR  (I )  ,CU  (I)  ,  PB  (I) 
INDX  ' 

IF 


NAME|J) 


NAM] 

DATA 


(J) 


15 


TO  APL  DATA  FILE 


(  (i  =j 

(  .GT.-99)  GO  TO  10 

CONTINUE 
WHITE  WEARMETAL 
WHITE  (10.*) 

WHITE  (9,*) 

WRITE  (8  ,  *) 

WRITE  (7  ,  *) 

WHITE  J11,*) 

GENERATE  HISTOGRAMS  OF  WEARMETAL  LEVELS 
CALL  HIST G?  (FE,I,15) 

WRITE  (6,200) 

CALL  HISTGP  (AL,I,15) 

WRITE  (6,210) 

CALL  HISTGP  (CR,I,15) 

WRITE  (6 ,220) 

CALL  HIST GP  (CU,I,15) 


LEVELS 
PB 
CU 
CR 
AL 
FE 


WRITE  (6,230) 

/■  i  r  r  ri"r  r»  m  «« 


CALL  Hi S ^ GP 
WRITE 


100 

110 

200 

210 

220 

230 

240 

245 

C250 

C270 


WHITE 

WRITE 

WRITE 

WHITE 

WRITE 

STOP 

FORMAT 

FORMAT 

FORMAT 

FORMAT 

FORMAT 

FORMAT 

FORMAT 

FORMAT 

FORMAT 

FORMAT 

END 


6,24  0) 

.g,* 

i6**) 


(EE, I, 15) 


FE 

AL 

CB 

CU 

PB 


4,0L 


6  FID .  1) 

IX,  1  HISTOGRAM 
IX,  »  HISTOGRAM 
IX,  •  HISTOGRAM 
IX,’  HISTOGRAM 
IX,  «  HISTOGRAM 
1  X,5F  10.5) 

IX,  6F  1C.  1) 
Ix'AIQ) 


OF 

OF 

OF 

OF 

OF 


PPM 

PPM 

PPM 

PPM 

PPM 


LEVELS 

LEVELS 

LEVELS 

LEVELS 

LEVELS 


FOR 

FOH 

FOR 

FOE 

FOR 


IRON'  ) 
ALUMINUM • 
CHROMIUM’ 
COPPER’) 
LEAD’) 


C  *****  ERCGRAM  ANALYSIS******* 

C  PROGRAM  ANALYZES  THE  10  HOUR  TRENDS  IN  WEARMETAL 
C  RICHARE  I.  BAUER 

REAL  DELHR(14QQ),  DEL  FE  (1400).DSLAL(  1400)  , 

♦  DELCR  <14  001  ,  DELCU  (14  00)  .DELPB  (14  00)  . 

♦  TRND  AL ( 1  4 00)  ,  TRNDCR  (  1400)  ,  TRNDPB  ( 1400)  , 

♦  HR  (1400)  ,FE (1400)  ,AL  (14  00)  ,CR  ' 

♦  FEJ1400)  .CUO'iOO)  ,  TR  NDCU  ( 1  400 
INTEGER  PRE 


TRNDPB  ( 1400)  , 
rCR (1400)  , 

400) , TRNDrE (1400) 


'Sfil 


DIMENSION  NAME  (75) 

1=1 
K=  1 

DO  50  J  =  1  ,75 
READ  (5,110)  NAME  (J) 
WRITE  (6,111)  NAME (j) 
CONTINUE 

READ  (5,100)  HR  (I) 
WRIT  E  (6,200)  HR(I 

*PB  (I) 

IF  (HR (I) .EQ.  (0  .  )  ) 
CONTINUE 

1=1+1 

CONTINUE 

READ  (5,100)  HR  (I) 
WRITE  (6,200)  HR<2 

*PB  (I) 

IF  (HR (I)  .  EQ.-99.) 
ERE  =  1-1 

IF  (HR  (I)  .EQ.  (0  .)  ) 


»FE  (I) 
),FE(i; 


AL(I)  ,CR(I)  ,CU  (I)  ,  PB  (I) 
,  AL  (I)  ,CS  (I)  ,CU  (I)  , 


5.100) 

(6 , 200] 


-99.)  GO  TO  49 


(HR (I)  .EQ.  (0  .))  GO  TO  7 
(HR  (I)  .LE.HR  (PRE)  )  GO  TO 
DELHH  (K)  =H  R  (I )  -HR  (PRE) 


WRITE 
*  EELCU 


DELI  E 
DEL  AL 
DEL  05 
DELCU 
DELEE  _ 
(6,*)  DELdR(K) 
(K)  ,  DELPEJK) 

GO  1(3  15 
CONTINUE 
DELHR(K)  =H 
DELFE(K)  =F 
DELAL(K)  =  A 


i  -FE (PR E 
l  -AL  (PRE 
-CR  (PRE 
l  -CD  (PRE 
- PB(PREl 


[)  ,  DEL  AL  (K)  ,  DELCR  (K)  , 


DELCR  (K 
DELC  U  (K 
DELPB  (K 
WRITS  (6 
*DELCU  (K)  ,  DELP 
CONTINUE 
TRNDEE(K) =(DE 
TEND AL(K  =  (DE 
TRNDCR (K  *i  DE 
TRNDOU (K  =  ( DE 


V  I  ' 

!I§U) 


DELHH  (K)  ,  DSLFE  (K)  ,  DEL  AL  (K)  ,  DELCR  (K)  , 


TRNDFE(K) =  < 


DEIFE 
DELALi 
DELCRl 
DELCU  i 
DELPB  i 
TRNDI 


WRITE  (6,*)  TRNDF  E  (K 
*TRNDEE(K) 

K=K+  1 
GO  TO  6 
CONTINUE 
1=1-1 
CONTINUE 

CALL  HISTGP  (TRNDFE, 
WRITE  (6,210) 

CALL  HISTGP  (TBNDAL, 
WRITE  (6 ,220) 

CALL  HISTGP  (TRNDCR, 
WRITE  (6 ,230) 

CALL  HISTGP  (TRNDCU, 


/DELHR 
/DELHH i 
/DELHR 
/DELHR 
/DELHR i 


*10. 
*  10. 
*10. 
*  10. 
*10. 


E(K)  , T RNB AL  (K)  , TRNDCR  (K)  ,  TRNDCU  (K)  , 


WRITE  (6,240) 


(TRNDCU, 


K,0) 

K.O) 

K,0) 

K,0) 


onoonnnnnn 


250 


FORMAT  (1X,6F10.5) 

FCRHAT  IX,  A10) 

FORMAT  1X,A10[ 

FORMAT  1  X,  6  F  ID.  2) 

FORMAT  1 X, *  HISTOGRAMS  OF  PIEC2WIS  2  TRENDS  FOR 
♦IRON*  ) 

FORMAT  (IX, 'HISTOGRAMS  OF  PIECEWISE  TRENDS  FOR 

♦  A1CMINUM '  ) 

FORMAT  (IX,  «  HISTOGRAMS  OF  PIECEWISE  TRENDS  FOR 
♦CCPPER* ) 

FORMAT  (IX, '  HISTOGRAMS  OF  PIECEWISE  TRENDS  FOR 

♦  CHROMIUM') 

PORMAT  (IX,' HISTOGRAMS  OF  PIECEWISE  TRENDS  FOR 
♦LEAD') 

END 


50 


c  *****  EHCGRAM  an algen  ******* 

C  PROGRAM  TO  GENERATE  SIMULATED  DATA  FOR  USE  IN  ANALYSIS 
C  RICHARE  F .  BAUER 

BEAL  HR  (600)  .IRON  (60  0)  ,  ERROR  l 240  0)  ,  DELHR  (599)  , 

♦  DELI RN (59  9)  , TRDIR  N  (5  99f  .IRONT  (24  00)  , 

♦  DELIB1J2396)  -TBDIRT  (  2396) 

INTEGER  PRE,PEEI,PREJ,PRr!K,PREL 
SEED=10. 

CALL  LNORMISEED, ERROR, 2400, 2,0) 

DC  5  1=1,550 

HR(I)=1Q.*  (FLO  A  T  (I)  ) 

IRON  (11=  (  (.3648)  *HR  (I))  +ERROR  (I) 


IRON  (11  = 

J=  I  ♦  60  Q 
K=  I  ♦  1 2  0  0 
L=I+  1800 

IRON1  (I)  =  (.  364  8)  *HR  (I)  ♦ ERROR  (I) 
IRON  1  (J  =  (.  4277  *HRiI  +ERROR  J 
IRON1/K  =11.5141  *HR  I  ♦  ERROR  (K 
IRON1(L  =i(.  330  6  *HR  I  +  ERROR  (L) 
CONTINUE 
EO=IRON (I) 

EOI=IRON1  (I) 

B  C J=IRON 1  i  J 
EO  K=  IRO  N  1  K 
BOL=IRON  1  iL 
DC  7  1=551,600 
JJ  =1-500 

HR(I)  =  1Q.*  (FLO  A  T  (J  J)  ) 

J=I*600 
K=  I  ♦  1 2  0  0 
L=I+  1800 

IRON  (I)  =EQ+ (  (.7  296)  *HR(I)  ) +ERROR  (I) 
IRO  N  i  (I)  =E0I+((.7296)  *HR(I)  )  +  ERR(JR| 


,  *  (FLOAT  (JJ)  ) 


IRO  N  1  ( J  = 
IRON1/K  = 
IRO  N  1  (L  = 
CONTINUE 
DC  1C  JJ  =2,  600 
I=JJ 

J=J  J+60  0 
K=JJ+1200 
L=JJ  +  1  800 
PRE= JJ- 1 


=E0l  +  ( ( . 7296)  *HR(I|  ) 
=  E0J  ♦  i'  f .9554)*HR  (in 
=E0K+((  1.0282)*HR(I) 
=E0L  +  i'  ( .66  12)  *HR  (I) 


♦ERROR  (I) 
♦ERRORjJ) 

)  +  ERROR (K) 
♦  ERROR  (L) 


DELHR(JJ)=HR  (JJ)  -HR  (PfiE) 

DELIRN  (JJ)=IRON  (JJ)  -IRON  (PRE) 

TRDIRN  (JJ)  =  (DELIRN  (JJ) /DELHR  (JJ)  )  *10. 


C200 


PRE  1=1-  1 
PR  E  J=J-  1 
PREK=K-  1 
PR  E  L  =L  -  1 
DEL IR1 (I  =IRC 
DELIR1  J  i  =IRC 
DSLIR1  (K  =IRC 
DELIR1  (L  =IRC 
TRDIR1  iT  =  (Di 
TRDIR1  i'J'i  =  (DF 
TRDIR1  i'  K  =(DI 
TRDIR1  iL  =(DI 
WRITE  i  6  ,200) 
CONTINUE 

CALL  HISTGP  (TRDIE 
WRITE  (4 ,  *1  TRDIRK 
CALL  HISTGP  (TRDIE 
STOP 

FORMAT  (1 X, IF  10.5) 
END 


=IRON 1  (I 
=IR0N1  i'J 
=IRON  1  (K 
=IR0N1  (L 
=  (DELIR1 
=  (DELI  R1 
= (DELI R1 
=  (DELIR1 
200)  TRD] 


-IRON  1  (P 
-IRON1 (P 
-IRON  1  (P 
-IRON  1  (P 
I)  /DELHR 
J  /DELHR 
’K'/DELHR 
L  /DELHR 
RN (JJ) 


PR  El) 
PR  EJ 
PEEK 
PEEL 
R  (JJ 
R  ( JJ 


*10. 

*10. 

*10. 

*10. 


(TRDIRN,  PRE, 15) 
TRDIRN 

(TBDIR1  ,  2396,15) 
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